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Integrated analyses identify the involvement of 
microRNA-26a in epithelial-mesenchymal transition 
during idiopathic pulmonary fibrosis 

H Liang 1 ' 3 , Y Gu 2 ' 3 , T Li 1 , Y Zhang 1 , L Huangfu 1 , M Hu 1 , D Zhao 1 , Y Chen 1 , S Liu 1 , Y Dong 1 , X Li 1 , Y Lu 1 , B Yang*' 1 and H Shan*' 1 

Idiopathic Pulmonary Fibrosis (IPF) is a chronic, progressive, and highly lethal fibrotic lung disease with poor treatment and 
unknown etiology. Emerging evidence suggests that epithelial-mesenchymal transition (EMT) has an important role in repair 
and scar formation following epithelial injury during pulmonary fibrosis. Although some miRNAs have been shown to be 
dysregulated in the pathophysiological processes of IPF, limited studies have payed attention on the participation of miRNAs in 
EMT in lung fibrosis. In our study, we identified and constructed a regulation network of differentially expressed IPF miRNAs and 
EMT genes. Additionally, we found the downregulation of miR-26a in mice with experimental pulmonary fibrosis. Further studies 
showed that miR-26a regulated HMGA2, which is a key factor in the process of EMT and had the maximum number of regulating 
miRNAs in the regulation network. More importantly, inhibition of miR-26a resulted in lung epithelial cells transforming into 
myofibroblasts in vitro and in vivo, whereas forced expression of miR-26a alleviated TGF-/M- and BLM-induced EMT in A549 cells 
and in mice, respectively. Taken together, our study deciphered the essential role of miR-26a in the pathogenesis of EMT in 
pulmonary fibrosis, and suggests that miR-26a may be a potential therapeutic target for IPF. 
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Idiopathic Pulmonary Fibrosis (IPF) is a chronic, irreversible 
and lethal lung disorder with unknown etiology. 1 It has been 
proposed that IPF likely results from the aberrant activation of 
alveolar epithelial cells after injury along with the accumulation 
of inflammation cells, which provoke the excessive prolifera- 
tion of fibroblasts and exaggerated deposition of extracellular 
matrix (ECM) in the lung parenchyma, and finally leading to 
the destruction of lung parenchyma. 2,3 Although significant 
progress has been made in recent years in elucidating key 
aspects of the pathobiology of IPF, the disorder lacks effective 
treatment strategies. 4 

Although the molecular mechanisms underlying IPF are 
largely unknown, many scientists believe that the fibroblasts 
foci formation and excessive ECM production are direct 
causes of IPF. 5,6 The key mesenchymal feature of patholo- 
gical fibrosis is the increased numbers of transdifferentiated 
fibroblasts, namely, myofibroblasts. Myofibroblasts share 
features with both fibroblasts and smooth muscle cells; they 
overexpress a-smooth muscle actin (a-SMA) and are most 
likely responsible for the enhanced synthesis of abnormal 
matrix observed in pulmonary fibrosis. 7 Myofibroblasts have 
multiple origins in the lung, such as fibroblasts, epithelial cells 



and bone marrow-derived circulating fibrocytes. Fibroblasts 
and epithelial cells are the primary origins of myofibroblasts in 
IPF. 8,9 According to lineage tracing in murine models of lung 
fibrosis, it is estimated that up to one-third of fibroblasts may 
be of epithelial origin. 10 ' 11 However, the role of epithelial- 
mesenchymal transition (EMT) in lung fibrosis is just begin- 
ning to be understood. Therefore, determining the mechanism 
of EMT involved in the pathogenesis of lung fibrosis is likely to 
suggest important new directions for the treatment of IPF. 

MicroRNAs (miRNAs) are a class of 19-22nt non-coding 
small RNAs that control the expression of a large number of 
genes by binding to the 3'UTR of targets and blocking 
translation or by causing degradation of target mRNA. 12 
miRNAs have been implicated in the initiation and progression 
of pathological processes, including cancer and fibrosis in the 
heart, liver, kidney and others. 13-15 Some papers have 
reported that miRNAs are involved in cancer by regulating 
the EMT. 16,17 However, limited studies have payed attention 
on the participation of miRNAs on EMT in lung fibrosis. 

In our study, we focused on investigating the impact and 
mechanism of miR-26a on EMT involved in IPF and the 
potential effect on improving treatment for IPF. Our results 
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showed that miR-26a was significantly down regulated in IPF, 
and the EMT pathway was significantly enriched with the 
differentially expressed IPF genes. Moreover, inhibition of 
miR-26a caused apparent EMT in A549 cells and in mice, 
whereas overexpression of miR-26a retarded Transforming 
growth factor (TGF)-/?1- and bleomycin (BLM)-induced EMT 
in vitro and in vivo, respectively. These findings indicate that 
miR-26a is a potential EMT inhibitor and may be considered a 
new therapeutic tool for IPF. 

Results 

EMT pathway is dysregulated in IPF. We detected 
significantly differentially expressed genes in the mRNA 
expression profile of IPF compared with the control group 
and performed functional enrichment using biological pro- 
cess terms in the Gene Ontology database. Our results 
showed that the EMT pathway (GO:0001837) was signifi- 
cantly enriched with the differentially expressed IPF genes 
(P=9.01E-4). The EMT pathway in the Gene Ontology 
database contains 79 genes, and 50 of them were detected 
in the microarray. We found that 35 out of the 50 detected 
genes were significantly differentially expressed in IPF. Two 
distinctive expression clusters of 35 EMT pathway genes 
were obtained by hierarchical clustering analysis (Figure 1). 
These results demonstrate that the EMT pathway plays 
important roles in the progression of IPF. 

To explore the role of miRNAs on regulating EMT in the 
process of IPF, we identified the differentially expressed 
miRNAs from the IPF miRNA expression profile compared 
with the control group and constructed a regulation network of 
differentially expressed IPF miRNAs and EMT genes 
(Figure 2a). The network contains all EMT pathway genes, 
including the differentially expressed genes, the non-differen- 
tially expressed genes and the genes not detected in the 
mRNA expression profile. We found that high-mobility group 



protein A2 (HMGA2) was regulated by 33 miRNAs that were 
differentially expressed in IPF (Figure 2b), which was the 
maximum number of regulating miRNAs, although HMGA2 
was not detected in the mRNA expression profile. This result 
indicates that HMGA2 may be dysregulated by miRNAs and 
may participate in the progression of EMT, which further 
induces IPF. Moreover, miR-26a showed a significant down- 
regulation in the IPF miRNA expression profile (P=7.17E- 
13), and miR-26a regulated many EMT-related genes 
(Figure 2c), including HMGA2. 

miR-26a is downregulated in mice with experimental 
pulmonary fibrosis. To investigate whether miR-26a parti- 
cipates in the development of pulmonary fibrosis, we 
constructed a mouse model of pulmonary fibrosis by 
administering intratracheal BLM, a well-studied pulmonary 
fibrosis model. 18 After 4 weeks, there was a significant 
increase in the deposition and content of collagen in the mice 
treated with BLM compared with the mice treated with saline 
(Figures 3a and b). Additionally, the mRNA levels of 
collagenlal and collagen3a1 were markedly upregulated in 
BLM-treated mice (Figure 3c). Real-time PCR analysis found 
that the expression of miR-26a was significantly decreased in 
the lungs of mice following administration of BLM (Figure 3d). 

It has been reported that TGF-/?1 -induced EMT is the major 
source of myofibroblasts, which has a substantial role in IPF 
by secreting ECM components to promote lung fibrosis. In our 
study, we found that both TGF-/?1 and p-Smad3 were 
upregulated at the protein level in the mice receiving BLM. 
In addition, the expression of HMGA2, a key factor in EMT, 
was strikingly elevated in the BLM models (Figure 3e). 
Furthermore, the expression of E-cadherin was markedly 
decreased, whereas the expression of a-SMA was signifi- 
cantly increased in the lungs of mice with BLM-induced 
pulmonary fibrosis. These results indicated that EMT 
occurred in our model of experimental pulmonary fibrosis. 
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Figure 1 Differentially expressed EMT pathway genes between IPF and control group. Heat map represents the differentially expressed genes (adjusted P-value < 0.01 ) 
of EMT pathways in IPF group compared with the control group 
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Figure 2 Regulation network of IPF miRNAs and EMT-related genes, (a) The regulation network of differentially expressed miRNA in IPF and EMT pathway genes, 
(b) The subnetwork of regulation relationship among the miRNAs and HMGA2, which has the maximum number of regulating miRNAs in the regulation network, (c) The 
network of regulation relationship among miR-26a and its targets. The size of each node corresponds to the number of genes or miRNAs connected to the node 



HMGA2 is a direct target of miR-26a in A549 cells. Based 
on the above findings, we proposed that miR-26a is involved 
in the process of IPF, and the negative correlation between 
miR-26a and EMT may participate in the progression of 
pulmonary fibrosis. Alignment of miR-26a with the HMGA2 
3'-UTR sequence revealed one potential conserved seed 
site, which appeared in several databases (TargetScan, 
miRDB, miRanda and PicTar; Figure 4a). As depicted in 
Figure 4b, transient introduction of miR-26a into A549 cells 
led to a significant decrease in HMGA2 at the protein level. 
Moreover, this inhibition was reversed by co-transfection with 
AMO-26a, and the negative control (NC) had no significant 
effect. The relationship between miR-26a and HMGA2 was 
further verified in Figure 4c. Co-transfection of miR-26a with 
the luciferase reporter carrying a portion of the human 
HMGA2 3'-UTR caused a significant decrease in luciferase 
activities compared with transfection with luciferase vector 
alone. The reduction of luciferase activities was efficiently 
reversed by AMO-26a (10nmol/l). Furthermore, we used an 
HMGA2 3'-UTR construct that was mutated at the predicted 
miR-26a site. miR-26a failed to diminish the mutated 
luciferase activities, whereas the mutated miR-26a was able 
to do so. Finally, real-time PCR showed that miR-26a had no 
effect on HMGA2 mRNA levels (Figure 4d) demonstrating 
that miR-26a directly regulates HMGA2 in a post-transcrip- 
tional manner. 

Loss of function of miR-26a induces EMT in vitro and 
in vivo. HMGA2 is a key positive regulator in the process of 
EMT and induces pulmonary fibrosis in mice. 19 The direct 
regulation of miR-26a on HMGA2 demonstrates that the 
downregulation of miR-26a may contribute to IPF by inducing 



EMT. Moreover, miR-26a expression was significantly 
suppressed by treatment of A549 cells with recombinant 
TGF-/?1, a classic EMT model in epithelial cells (Figure 5a). 
To determine the potential role of miR-26a in the course of 
pulmonary fibrosis, we transfected AMO-26a into A549 cells. 
As shown in Figure 5b, AMO-26a effectively inhibited the 
expression of miR-26a and caused the elevation of collagen 
content (Figure 5c). Furthermore, the reduction of miR-26a 
induced an obvious EMT phenotype in A549 cells (Figure 5d) 
together with clear alteration of certain related proteins, 
including E-cadherin, vimentin and a-SMA (Figure 5e). 
Furthermore, we observed apparent EMT and lung fibrosis 
in the mice treated with an intratracheal injection of 
antagomiR-26a (Figure 6), which were observations similar 
to those detected in vitro. These results suggest that loss of 
function of miR-26a could facilitate lung epithelial cells to 
transform into myofibroblasts and induce pulmonary fibrosis 
in mice. 

Forced expression of miR-26a diminishes EMT and 
fibrogenesis in A549 cells and in mice with experimental 
pulmonary fibrosis. To evaluate the effects of miR-26a on 
EMT and lung fibrosis, we examined whether forced 
expression of miR-26a could mitigate the phenotype of 
EMT induced by TGF-^ in A549 cells. We transfected miR- 
26a into cultured A549 cells to assess the effects of miR-26a 
on collagen production. As shown in Figure 7a, miR-26a 
abolished TGF-/?1 -induced secretion of collagen in A549 
cells, while the NC group failed to do so. Additionally, 
immunofluorescence and Western blot assays showed that 
miR-26a attenuated TGF-^1 -induced EMT in A549 cells 
(Figures 7b and c). 
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Figure 3 miR-26a is downregulated in lungs from BLM -treated mice, (a) HE and Masson Trichrome staining of mouse lung sections shows interstitial fibrosis with collagen 
deposition 28 days after injection of BLM (magnification, 200 x ). (b) Increased collagen content in the lungs of mice treated with BLM compared with saline group, 
(c) Increased collagen I and collagen III mRNA expression in the BLM-treated mice compared with control animals, measured by real-time PCR. (d) Real-time PCR analysis 
of decreased expression of miR-26a in the lungs of mice with experimental pulmonary fibrosis, (e) Western blot analysis of proteins in the TGF-/? signaling pathway and 
EMT-related genes in BLM-treated mice as compared with control animals. Mean ± S.E.M.; n = 5 mice in each group; *P<0.05, **P<0.01 versus saline group 



To examine whether miR-26a has protective roles in the 
progression of lung fibrosis, we administrated an intratracheal 
injection of agomiR-26a 7 days after instillation of BLM and 
analyzed the extent of EMT and fibrosis 28 days after 
administration of BLM. As shown in Figures 7d-f, agomiR- 
26a significantly suppressed collagen content as well as the 
occurrence of EMT in BLM-treated lungs. 

Discussion 

In the present study, we took bioinformatic methods to identify 
miRNAs potentially involved in EMT under condition of 



pulmonary fibrosis. Combined with BLM-induced model of 
pulmonary fibrosis in mice, we confirmed the downregulation 
of miR-26a in fibrotic lungs. This reduction correlated with 
increased expression of mesenchymal cells and EMT-related 
genes. Further studies showed that HMGA2, a key positive 
regulatory factor in EMT, is one of the direct targets of miR- 
26a. Moreover, miR-26a could be suppressed by TGF-/?1 in 
A549 cells, and knockdown of miR-26a caused EMT in vitro 
and in vivo. Interestingly, forced expression of miR-26a 
diminished EMT and fibrogenesis in A549 cells and in mice 
with experimental pulmonary fibrosis. Taken together, our 
study supports a role for miR-26a in the pathogenesis of EMT 



Cell Death and Disease 



miR-26a regulates EMT in IPF 

H Liang et al 



a Mouse miR-26a_ 

3' ~- ucggauaggaccuaSugaacuu"]- 5' 

I I I llllllll 

1879 — CCAACCUUCAAAGAAACyyGAAI---1900 
Mouse HMGA2 3' -UTR (NM_01044l".2) 

Rat miR-26a_ 

3' - ucggauaggaccuaKugaacuu"!--- 5' 

II i llllllll 

1921 -- CCAACCUUCAAAGAAACJJUGAAJ--- 1942 
Rat HMGA2 3' -UTR (NM_032070.1) 

Human miR-26a 

3' - ucggauaggaccuaKugaacuu"]-- 5' 

ii iiiiiiii! 

1851 — CCAACCUUCAAAGAAACUUGAAI-1872 
Human HMGA2 3' -UTR (NM_003483.4) 



HMGA2 1 



GAPDH 



18kDa 



37 kDa 



* 

!## 
■iii 



Luciferase Assay 



oiiiiili 



«B 0.5] 



0.0 1 



HMGA2 mRNA 



buii 



Figure 4 miR-26a post-transcriptionally regulates HMGA2. (a) Sequence alignment showing miR-26a/HMGA2 complementarity for mouse, rat and human genes. The 
matched base pairs are outlined by dashed red rectangles. The Genbank accession numbers of the genes are indicated in the brackets, and the positions of the target sites are 
numbered, (b) Compared with control, transfection of miR-26a results in a significant decrease of HMGA2. Co-application of miR-26a with AMO-26a alleviated the reduction of 
HMGA2, whereas NC showed no effects, (c) Compared with control, transfection of miR-26a with the luciferase reporter gene vector containing the wild-type 3' UTR of 
HMGA2 results in a significant decrease of luciferase activity. Co-application of miR-26a with AMO-26a alleviated the reduction of luciferase activity, whereas NC showed no 
effects, (d) Real-time PCR shows that miR-26a had no effects on HMGA2 mRNA level. AMO/miR-26a inhibitor, AMO-26a; NC/negative control. Mean ± S.E.M.; n = 4 
represents three independent experiment under each condition; *P<0.05, **P<0.01 versus control; # P<0.05, ## P<0.01 versus miR-26a, & P<0.05 versus mis-vector 



in pulmonary fibrosis, and suggests that miR-26a may be a 
potential therapeutic target for IPF. 

TGF-/71 has been shown to have key roles in pulmonary 
fibrosis, not only through stimulating proliferation of fibro- 
blasts, but also through induction of EMT in alveolar epithelial 
cells by activating Smad- or non-Smad signaling pathways. 20 
The EMT process is characterized by the loss of epithelial cell 
phenotypes, such as diminished expression of E-cadherin 
and ZO-1, and the acquisition of mesenchymal cell pheno- 
types, such as enhanced expression of vimentin, a-SMA and 
N-cadherin. EMT converts adherent epithelial cells to motile 
mesenchymal cells with enhanced migratory capacity, inva- 
siveness, elevated resistance to apoptosis and greatly 
increased production of ECM components. 8 Emerging evi- 
dence suggests that EMT participates in embryonic develop- 
ment, tissue repair and pathological processes, notably tumor 
invasiveness and metastasis. 21 

Although large number of studies have confirmed the 
mechanism of EMT in different cancer types, little is known 
regarding the role of EMT in IPF. Current studies show that the 
role of EMT in the formation of activated myofibroblasts during 
the process of pulmonary fibrosis is still controversial. 22 ' 23 
Yamada et al. 24 did not find histological evidence for EMT in 
tissue samples from patients diagnosed with IPF and non- 
specific interstitial pneumonia using dual immunohistochem- 
istry. Roch et al. 25 did not find evidence at the cellular or 
molecular level for the EMT of labeled cells into myofibro- 
blasts by analyzing normal and fibrotic mouse and human 
lungs. On the contrary, an increasing number of studies have 
identified a major role of EMT in lung fibroblasts accumulation 



in IPF. A recent study by Chen YL et al 28 reported that 
sorafenib ameliorates BLM-induced pulmonary fibrosis along 
with the inhibition of EMT and fibroblast activation. Sonnylal S 
et al. 29 found that overexpression of connective tissue growth 
factor (CTGF) in epithelial cells causes EMT-like morpholo- 
gical changes and expression of a-SMA, which in turn 
contribute to lung fibrosis. Balli D et al. found that foxml is 
required for radiation-induced pulmonary fibrosis by enhan- 
cing expression of the genes critical for lung inflammation and 
EMT. 30 Our results are consistent with these abovementioned 
findings that the EMT is likely an important mechanism for 
IPF. Our results presented here strongly suggest for the first 
time the role of EMT in IPF by large-scale approach in 
numerous IPF samples. Additionally, the regulation network 
constructed by the differentially expressed miRNAs and EMT- 
related genes in IPF should also aid better understanding of 
the mechanisms for IPF. 

To examine the role of miRNAs on EMT, many scientists 
have studied the effects of miRNAs on the direct transcription 
factors related to EMT. Chang et al? A found that miR-200c 
has a role in regulating both EMT and EMT-associated stem 
cell properties in breast cancer through directly regulating 
ZEB1. Yang et al. 32 confirmed that miR-506 is a potent EMT 
inhibitor and therapeutic tool in serious ovarian cancer by 
targeting Snail2. Previous studies have confirmed that 
HMGA2 is a critical factor required for TGF-/71 -induced EMT 
in mammary epithelial cells. 33 In mammary epithelial cells, 
TGF-/71, via the Smad pathway, induces the expression of 
HMGA2, which binds to the SnaiM promoter and induces 
SnaiM expression, E-cadherin repression and the overall EMT 
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Figure 5 Inhibition of miR-26a causes EMT in cultured A549 cells, (a) TGF-/71 inhibits the expression of miR-26a in cultured A549 cells, (b) Real-time PCR analysis of 
miR-26a expression in A549 cells transfection of AMO-26a as compared with control, (c) Increased collagen content in A549 cells transfected with AMO-26a. 
Immunofluorescence (d) and Western blot (e) examine the effect of AMO-26a on EMT. n = 4; *P<0.05, **P<0.01 versus control 



phenotype. 34,35 However, few studies have focused on the 
role of HMGA2 in the process of IPF. In our study, we found 
that HMGA2 had the maximum number of regulating miRNAs 
in IPF-related miRNA-EMT gene regulation network. There- 
fore, we supposed that influencing HMGA2 instead of its 
downstream effectors may be a more effective way to control 
the process of EMT. Our study revealed that HMGA2 
upregulation by miR-26a inhibition is a new mechanism for 
the EMT involved in pulmonary fibrosis. 

The screening and study of miRNA involvement in 
pulmonary fibrosis is still in the early stages. Two recently 
published studies have confirmed the effect of dysregulated 
expression of let-7d and miR-200 on EMT in pulmonary 
fibrosis. 3,36 Pandit KV et al. found that let-7d was significantly 



reduced in human IPF tissues, and the inhibition of let-7 in 
epithelial cell lines in vitro and in vivo resulted in increased 
expression of mesenchymal-specific genes and pulmonary 
fibrosis. 3 Studies from Yang et al. 36 confirmed the down- 
regulation of the miR-200 family in the lungs of mice with 
experimental lung fibrosis and patients with IPF. 

Accumulating evidence has confirmed that miR-26a is 
involved in the pathological processes of many diseases, 
such as cardiac fibrosis, atrial fibrillation, angiogenesis and 
others. 37-39 However, the role of miR-26a on IPF, espe- 
cially on EMT, remains largely unknown. In this study, we 
found dysregulated expression of miR-26a in IPF. In 
addition, we showed that miR-26a could regulate the 
expression of HMGA2 and attenuate EMT, which suggests 
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Figure 6 Inhibition of miR-26a in the lung promotes EMT in mice, (a) Real-time PCR analysis of miR-26a expression in mice treated with antagomiR-26a compared with 
the saline group, (b) Increased collagen content in mice treated with antagomiR-26a. Immunofluorescence (c) and Western blot (d) assays were used to assess the EMT 
phenotype in antagomiR-26a-treated mice compared with control animals. n = 4; **P<0.01 versus saline 



that miR-26a could serve as a mediator to prevent fibrotic 
lung diseases. 

Numerous studies have indicated that miR-26a has a 
protective role in the process of cardiac fibrosis by regulating 
the activation of cardiac fibroblasts. 40,41 Therefore, we assume 
that miR-26a also has an important role in modulating the lung 
fibroblast phenotype based on the following: (1) Similar 
physiological characteristics and pro-fibrotic pathways exist in 
lung and cardiac fibroblasts. (2) Our study found that HMGA2 
was a direct target of miR-26a in A549 cells. Huleihel L et al. 42 
confirmed that let-7d significantly attenuated the differentiation 
of fibroblasts into pathological myofibroblasts by targeting 
HMGA2. (3) The data obtained from databases (TargetScan, 
miRDB, miRanda and PicTar) showed that miR-26a has many 



potential targets involved in the regulation of fibroblast 
phenotype and function, such as CTGF, fibroblast growth 
factor receptor, hepatocyte growth factor and so forth. IPF is a 
complex pathological process characterized by excessive 
fibroblast proliferation and injury of alveolar epithelial cells. 
Hence, it is undeniable that the anti-fibrotic effect of miR-26a in 
IPF is a combined effect mediated by a variety of mechanisms. 
Further studies are needed to reveal other potential mechan- 
isms involved in the anti-fibrotic effects of miR-26a and assess 
the combined effects of miR-26a in IPF. 

In addition to revealing the role and mechanism of miR-26a 
in EMT, our study explored the potential therapeutic value of 
miR-26a for IPF. Although many miRNAs have been shown to 
have considerable roles in regulating numerous diseases, the 
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Figure 7 Enhancement of miR-26a expression attenuates EMT in vitro and in vivo, (a) miR-26a suppresses the TGF-/71 -induced increase in collagen content in A549 
cells. Western blot (b) and immunofluorescence (c) assays found that miR-26a suppresses TGF-/?1 -induced EMT in A549 cells, (d) Administration of agomiR-26a attenuates 
collagen content in bleomycin-treated mice, whereas agomiR-NC has no effect. Western blot (e) and immunofluorescence (f) assays shows that agomiR-26a attenuates EMT 
in the lungs of mice with pulmonary fibrosis. E-cadherin is stained in green, Vimentin is in red and the nuclei are in blue (magnification, 200 x ). n = 4; *P< 0.05 versus control, 
#P< 0.05 versus TGF-/21; **P< 0.01 versus saline, ## P< 0.01 versus BLM 
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complexity of delivering miRNA in vivo and pharmacokinetics 
have limited their application in human study. Future studies 
are needed to validate the clinical value of miR-26a for IPF 
treatment. 

Materials and Methods 

Microarray datasets and detection of differentially expressed 
miRNAs and genes. The microarray datasets were downloaded from the 
database of Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/). The 
mRNA expression profile consists of 119 IPF patient samples and 50 non- 
diseased control samples (GSE32537). The miRNA expression profile has 106 
IPF patient samples and 50 non-diseased control samples (GSE32538). 43 The 
quantile normalization and log2 transformation were applied to the raw data. The 
values of replicate probes corresponding to a miRNA or a gene were averaged to 
represent the expression value of the miRNA or gene. A two-sample f-test was 
performed to detect differentially expressed miRNAs and genes, respectively. The 
P values were adjusted by the Benjamini and Hochberg (BH) correction procedure 
to account for multiple tests with false discovery rate < 1%. 

miRNA targets databases and pathway information. Considering 
that miRNA targets predicted by multiple algorithms might be more reliable, we 
used the miRNA-target interactions of Homo sapiens appearing in at least two of 
the nine datasets (TargetScan, miRanda, PicTar, miRBase, DIANA-microT, PITA, 
RNAhybrid, miRTarBase and miRecords). The human genes annotated in the 
biological process of 'EMT' were downloaded from the Gene Ontology database 
(GO:0001837). The hypergeometric distribution model was used to test whether 
the EMT pathway was enriched with differentially expressed IPF genes. 
Hierarchical clustering of expression values was performed with R software using 
the metric of Euclidean distance and average linkage. 

Experimental pulmonary fibrosis model and treatment. Animals 
(8-week-old C57BL/6 mice; 20-30 g) in this work were used in accordance with the 
regulations of the Ethic Committees of the Harbin Medical University and 
conformed to the Guide for the Care and Use of Laboratory Animals published by 
the US National Institutes of Health (NIH Publication No. 85-23, revised 1996). 
The mice were anesthetized with pentobarbital (40mg/kg, i.v.) and randomly 
divided into two groups: saline and BLM. To induce pulmonary fibrosis, BLM was 
intratracheal ly injected at a dose of 1 .5 U/kg body weight on Day 0, and the mice 
were killed at 4 weeks after the BLM treatment. The control group received the 
same volume of sterile saline. All surgical procedures were performed under sterile 
conditions. The cholesterol-conjugated miR-26a mimics/antisense and NC 
(agomiR-26a/antagomiR-26a and agomir-NC/antagomiR-NC, respectively) were 
purchased from RiboBio Co., Ltd. (Guangzhou, China). Mice were treated either 
with agomiR-26a/antagomiR-26a (10mg/kg body weight in 50 saline) or 
agomiR-NC/antagomiR-NC by intratracheal injection. 

Masson's trichrome staining. The lungs of mice in the BLM and saline 
groups were quickly dissected and immersed in 4% neutral buffered formalin for 
24 h and stained with Masson's trichrome. Fibrosis tissue was quantified with 
Image Pro Plus 6.0 (MediCybernetics, Silver Spring, MA, USA). Three areas were 
analyzed in each slide, and each area was divided into 100 squares. The number 
of collagen points (blue stain) was scored as 1 (present) or 0 (absent). The results 
are shown as the percentage of area occupied by fibrosis to the total area. 

Collagen content assay. Total collagen content was measured using the 
Sircol Collagen Assay Kit (Biocolor, Belfast, Northern Ireland) according to the 
manufacturer's instruction. Each sample was incubated at 4°C with stirring for 
24 h. After centrifugation, 100 /A of each supernatant was assayed. One milliliter of 
Sircol dye reagent bound to collagen was added to each sample, and the solution 
was mixed at 4°C for 30min. After centrifugation, the pellet was suspended in 
700^1 of the alkali reagent included in the kit and read at 540 nm with a 
spectrophotometer. Collagen standard solution was utilized to construct a 
standard curve. Total collagen content was determined by assaying cell protein 
content. The amount of protein was colorimetrically assessed at 562 nm with a 
BCA protein assay kit (Beyotime, Nantong, China). 

Real-time RT-PCR. Total RNA samples were extracted from the lung tissues 
of mice or A549 cells using Trizol (Invitrogen, Carlsbad, CA, USA). As delineated 



in the work of Yang etai, the relative expression levels of mRNAs and miRNAs 
were quantified by the mirVana qRT-PCR miRNA Detection Kit (Applied 
Biosystems, Foster City, CA, USA), and real-time RT-PCR with TaqMan 
microRNA assays (Applied Biosystems), respectively. After the cycle reactions, 
the threshold cycle (Ct) was determined and relative mRNA and miRNA levels 
were calculated based on the Ct values and normalized to GAPDH or U6 levels in 
each sample. 

Cell culture and treatment. A549 cells were purchased from the Cell Bank 
of Chinese Academy of Sciences (Shanghai, China). The cells were grown in 25- 
cm 2 cell culture flasks with F-12K medium (Hyclone, Logan, UT, USA) containing 
10% fetal bovine serum (Hyclone), 2mM L-glutamate, 100 U/ml penicillin G, and 
100 U/ml streptomycin at 37 °C in 5% C0 2 , 95% air. The cells were maintained at 
sub-confluent densities and used for experiments when reaching passages 3 or 4. 
After starvation in serum-free medium for 24 h, A549 cells were administered 
recombinant human TGF-/71 (Sigma-Aldrich, St. Louis, MO, USA) for 48 h. 

Transfection procedures. For transfections, the cells were washed with 
serum-free medium once and then incubated with 4 ml of serum-free medium for 
4-6 h. The miRNA and Lipofectamine 2000 (Invitrogen) were separately mixed 
with 500 yul of Opti-MEM I Reduced Serum Medium (Gibco, Grand Island, NY, 
USA) for 5 min. Then, the two mixtures were combined and incubated at room 
temperature for 20 min. The Lipofectamine-miRNA mixture was added to the cells 
and incubated at 37 °C for 6h. Subsequently, 4 ml of fresh medium containing 
10% FBS was added to the flasks and the cells were maintained in the culture until 
the following experiments. 

Immunofluorescence staining. Cultured A549 cells on sterile glass cover 
slips were briefly washed with cold PBS 3 times and fixed with 4% 
paraformaldehyde for 15 min. Then, the cell membrane was penetrated by Triton 
X-100 for 1 h and blocked by normal goat serum for 1 h at 37 °C. The cells were 
incubated with anti-E-cadherin or anti-vimentin antibody (Abeam Inc., Cambridge, 
MA, USA; Cell Signaling Technology, Beverly, MA, USA) overnight at 4 °C and 
subsequently incubated with FITC-conjugated goat anti-mouse or goat anti-rabbit 
antibody for 1 h. The cells were washed with PBS, and the nuclei were stained 
with DAPI (Roche Molecular Biochemicals, Basel, Switzerland) for 5 min at room 
temperature. Immunofluorescence was analyzed under a fluorescence microscope 
(Nikon 80 i, Japan). 

Western blotting. For Western blot analysis, total protein samples were 
extracted from tissues or cells with the procedure as previously described. 45 Tissues 
or cells were lysed with RIPA lysis buffer (Beyotime). A total of 50 ^g of protein was 
fractionated on a 10% or 15% SDS-polyacrylamide gel. After electrophoretically 
transferring to a Pure Nitrocellulose Blotting membrane (Pall Life Sciences, Ann 
Arbor, Ml, USA), the blots were probed with primary antibodies, with GAPDH (anti- 
GAPDH antibody from Kangchen, Shanghai, China) as an internal control. Primary 
antibodies against TGF-/71, p-Smad3 and vimentin were rabbit polyclonal antibodies 
that were purchased from Cell Signaling Technology. Primary antibodies against 
E-cadherin or a-SMA were mouse polyclonal antibodies that were purchased from 
Abeam Inc. The immunoreactivity was detected using Odyssey Infrared Imaging 
System (Gene Company Limited, Hong Kong, China). The bands were quantified by 
measuring the band intensity for each group. 

Luciferase reporter assays. HMGA2 3'-UTR containing the conserved miR- 
26a binding sites was synthesized by Invitrogen and amplified by PCR. The PCR 
fragment was subcloned into the Sad and Hindi II sites downstream the luciferase 
gene in the pMIR-REPORT Luciferase vector (Applied Biosystems). The luciferase 
vector (100 ng) containing the 3' UTR was cotransfected with miR-26a mimics or 
AMO-26a into HEK-293 cells using Lipofectamine 2000 (Invitrogen). As an internal 
control, 10 ng of Renilla luciferase reporters were also included. After a 36-h 
transfection, the cells were collected, and the dual luciferase activities were measured 
using a luminometer (Promega, Madison, Wl, USA) according to the manufacturer's 
instructions. 

Statistical analyses. All data are presented as the mean ± S.E.M. One way 
analysis of variance ANOVA followed by Dunnett's test was used for multiple 
comparisons. A two-tailed value of P<0.05 was considered a significant 
difference. The data were analyzed using GraphPad Prism 5.0 (GraphPad Inc, 
San Diego, CA, USA) and SPSS 14.0 (SPSS Inc, Chicago, IL, USA). 
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